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The cross section of Xcj hadroproduction is calculated in 
the k± -factorization approach. We find a significant contri- 
bution of the Xci state due to non-applicability of the Landau- 
Yang theorem because of off-shell gluons. The results are in 
agreement with data and in contrast to the collinear factoriza- 
tion show a dominance of the color singlet part and a strong 
suppression of the color octet contribution. Our results could 
therefore lead to a solution of the longstanding controversy 
between the color singlet model and the color octet mecha- 
nism. 

The production of heavy quarkonia received a lot of at- 
tention from both theory and experiment in recent years. 
It is e.g. the most prominent signal in the search for the 
quark gluon plasma. Its usefuUness is, however, question- 
able as long as the charmonium production process is not 
understood. For a review we refer to Originally 
heavy quarkonium production was described in the color 
singlet model (GSM) QJ^. Calculations based on this 
model and standard collinear factorization show however 
disagreement with the experimental data. For example 
the next-to-leading order (NLO) QCD collinear results 
for direct J/^ hadroproduction underestimate the mea- 
sured cross section at Tevatron by a factor of « 50 (see 
fig. 4 in 1^ and Ref. 0). The proposed solution to this 
strong discrepancy is the so called color-octet-mechanism 
(COM) ||,^, according to which a color octet qq-pair 
which has been produced at short distances can evolve 
into a physical quarkonium state by radiating soft gluons. 
The COM introduces uncalculable non-perturbative pa- 
rameters, the color octet matrix elements, which have to 
be determined by a fit to the data ||l^,|ll| . The inclusion 
of the COM into NLO QCD collinear calculations leads 
in the case of hadroproduction to a reasonable agreement 
with experiment |l^Jll|. In these calculations the color 
octet contribution dominates. 

On the other hand up to now the COM suffers at 
least from two unsolved problems. When the, suppos- 
edly universal, color octet matrix elements are applied 
to electroproduction of heavy quarkonium the theoreti- 
cal predictions fail to describe the data . Furthermore 
the results of the COM for polarized heavy quarkonium 
hadroproduction seem to be incompatible with recent 
data from Tevatron Q . 

The longstanding discrepancy between the results 
based on the CSM together with collinear factorization 
and the experimental data shows up especially strongly in 



the /cj_-dependent cross sections from Tevatron |^]. Thus 
one can wonder if the collinear approximation, in which 
in NLO the only transverse momentum of the produced 
quarkonium comes from an additional final state gluon, 
is suitable at all. 

The aim of our paper is to clarify this question 
by a study of Xcj hadroproduction within the k±_ - 
factorization approach, which takes the nonvanishing 
transverse momenta of the colliding t-channel gluons into 
account. Generically this corresponds to taking into ac- 
count new regions of the phase space of the colliding 
gluons which is mandatory for the description of hard 
processes in the Regge region. 

More precisely we calculate the production of J/^f's 
originating from radiative Xcj decays. In a recent study 
iQ of open bb hadroproduction we found that k± - 
factorization gives far better results than NLO collinear 
QCD calculations and we expect a similar improvement 
for heavy quarkonium production. The main ingredients 
of our calculations in are the unintegrated gluon dis- 
tribution and the effective next-to-leading-logarithmic- 
approximation (NLLA) qq -BFKL production vertex 
which we use in this article as well. The projection of the 
heavy quark-antiquark pair onto the corresponding char- 
monium state is described in the standard way within 
the non-relativistic-quarkonium-model |^,^,|lO|,0 . 




FIG. 1. The basic diagram 

We study the production of Xcj whose lowest Fock 
state component is qqi^Pj). For J/'^ (a qq{^ Si) state) 
the LO production amplitude is zero. In order to get a 
nonzero (7g("^5i)-amplitude one has (in NLO in as) to 
emit an additional gluon. The amplitude for the produc- 
tion of a qq -pair plus a gluon within the BFKL approach 
would in our case require an effective three-particle pro- 
duction vertex which still has to be derived. In contrast 
the production of a Xci can be calculated in our approach 
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in LO because the Landau- Yang theorem which usually 
forbids the production of a Pi state is not valid for off- 
mass-shell gluons. 

We use the following definition of the light cone coor- 
dinates 

k+ = k°+k^, k- = k^^k^, ki_ = (0, fc\ fc^ 0) = (0, k,0). 

In the cm. frame the momenta of the scattering hadrons 
are given by 



The functions 6(^1,^2) and b'^{k2,ki) are illustrated 
in Fig.^ and their explicit form can be found in [Q. 
One important property of the charmonium production 
amplitude for on-mass-shell quark and antiquark states 
(^) , which is related to the gauge invariance of the whole 
approach, is its vanishing in the limit qi± —>■ (or q2±_ — » 
0). 



where the Mandelstam variable s is as usual the c.m.s. 
energy squared. The momenta of the t-channel gluons are 
qi and 52 (see Figjl]). The on-shell quark and antiquark 
(with mass m) have momentum fci respectively /c2 with 
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In the high energy (large s) regime we have 

P^=qi -q} ~ Qi-^ P^ = qi -q2 ~ -q2, 9i/2 ~ qi/2±^ 

where P = fci -I- fc2 is the momentum of the heavy quarko- 
nium with P^ = 4m^. The longitudinal momentum frac- 
tions of the gluons are xi = q^ /P^ , X2 = —q2 IP2 ■ 

The heavy quarkonium hadroproduction cross section 
in the /c^-factorization approach is [0, 

1 (d^P 



FIG. 2. The effective vertex 

The relation between the usual gluon distribution 
xg{x, q^) and the unintegrated gluon distribution J-{x, k) 
is given by 



xgix,q^) 



1^ 



e(q2-k2).F(x,k). 



(3) 



J-'(x, k) includes the evolution in x and k^ described 
by the BFKL and DGLAP equation. In the non- 
perturbative region of small k^ the unintegrated gluon 
distribution is not known, therefore we write (|^) accord- 
ing to llUlllllll as 



1 fd"^ P 

- W^Jtt'''^-'''^^-''^'^^- - - xgix,e) = ..9(.,qg) + /; ^e(q^ - k^)^(.,k) 



T{xi,qi±)- 



^ix2,q2±)- (1) 



The factor {N^ — 1)^ comes from the projection on color 
singlet in the t-channel. J-'{x,qj_) is the unintegrated 
gluon distribution. The heavy quarkonium production 
amplitude 4''x^'^^{xi^X2,qi±,q2±, P) is factorized (see be- 
low) in a hard part which describes the production of 
the qq pair and an amplitude describing the binding of 
this pair into a physical charmonium state. We choose 
the scale for asifJ-"^) in the amplitude ^p'^'^^ to be 
= respectively q^ = -qj^ 

The amplitude for the production of the charmonium 
state can be written as 



which introduces the a priori unknown initial scale qo and 
the initial gluon distribution xg{x, qg). Following pl| , p2t , 
we neglect the momentum dependence of the hard cross 
section in the soft region |q| < |qo|, so that 
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^T"" ^'Piqq^xc.i)*-^"'"'- 



(2) 



The qq production vertex vj/'^^ci (derived in j20| for mass- 
less QCD, appropriately generalized for massive quarks, 
has the form 



[^(<Zl±, g2±)e(q2 - q2) + 5(gi^, 0)e(q2 - q^)] e(q? - q^) 
+ [5(0, g2±)e(q^ - q2) + 5(0, 0)e(qg ~ q^)] e(q2 - q?), 

see also the discussion of this expression in . 

One important point is the proper choice of the unin- 
tegrated gluon distribution function. We use the results 
of Kwiecinski, Martin and Stasto They determined 
it using a combination of DGLAP and BFKL evolution 
equations. With the initial conditions 



qg = 1 GeV, xg(a;,q2) = 1.57(1 -x) 
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(4) 



- -.9' {t^'t^" h{ku k2) - t'H'' b^ik2, fci)) , 



where are the colour group generators in the fundamen- 
tal representation. The operator V{qq Xcj) projects 
the qq pair onto the charmonium bound state, see below. 



they obtained an execellent fit to P2(a;, Q^) data over a 
large range of x and Q^. 

In order to see the effect of off-shell gluons and the 
inapplicability of the Landau- Yang theorem as well as 
to perform calculations which do not require a fit to the 
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data we start with calculation of the color singlet part of 
the amplitude. This is most easily done by adapting the 
method of The projection of the hard amplitude 

onto the charmonium bound state is given by 



J + 7) = C(10 ^) the contribution from Xcj=o is neg- 
ligible. For the numerical computation we use the values 
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(L = 1, L,, ^ = 1, ^,1 J, J,) (3i, 3j|l) Tr {*^f ^7's=i,5. } , 

(5) 

where ?'l=i,l^ {Q = ki — is the momentum space wave 
function of the charmonium, and the projection operator 
Vs=i.s^ for a small relative momentum q = ki — k2 has 
the form 

The Clebsch-Gordan coefficient in color space is given by 
(3i,3j|l) = 5ji/\/Nc- Since P-waves vanish at the ori- 
gin, one has to expand the trace in in a Taylor series 
around q = 0. This yields an expression proportional to 



^(7"'i>L=i,L.(<?) - -^\/5e"(L,)7^'(0), 



(27r) 
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The cross section for J/^P production from radiative Xcj 
decays is then given by |10|,pT[ 



CTj/* from Xo = X! '^PiP2^XcjX ■ Br{XcJ ^ J/'^+l), 
J=0,l,2 

with the XcJ hadroproduction cross section op-^p^^^^jX 
(0). Because of the small branching ratio Br{xcj=o ^ 



mc = 1.48 GeV, |7^'(0)p = 0.075 GeV^ 
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M^ + Ps)/{y^lf^M^-P,)). To 



compare with data we multiply our cross sections with 
the braching ratio Br{J/'i> ji^ ^i^). 
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with the derivative of the P-wave radial wave function at 
the origin TV{Q) whose numerical values can be found in 
[^6|. For the individual Xcj=i and XcJ=2 amplitudes we 
use 

where we introduce the spin 1 and spin 2 polarization 
tensors ep{Jz) and e'^''{Jz) of the produced charmonium 
XcJ=i respectively XcJ=2- In the unpolarized case the 
squared amplitudes are further evaluated using 
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FIG. 3. The transverse momentum differential cross sec- 
tion in comparison to the data and a NLO QCD calculation 

The resulting Pi-dependent cross section for J/\I''s 
from radiative deacays of Xc's produced in pp-collisions 
is shown in Fig'^ together with the data from the CDF 
Collaboration ff and a NLO QCD coUinear result (see 
Fig. 7 in The individual contributions from Xci and 

Xc2 are shown in Fig.^. 



CS k,-factorization 
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FIG. 4. The individual contributions from Xci and Xc2 

The description of the data by the color singlet part 
alone is very satisfactory and becomes even better if the 
difference of the transverse momentum of J/^! (which 
is measured experimentally) and Xc (which enters our 
calculation) is taken into account. (Due to the radiative 
decay the transverse momentum of J/\l/ is typically larger 
by an amount of « 300 MeV than the corresponding 
Xc one which leads to a shift of the theoretical curve to 
the right.) The typical scale of the gluon off-shellness 
is given by the transverse momentum of the produced 
quarkonium. 
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We emphasize that the result has been obtained with- 
out fitting any of the parameters involved: The un- 
integrated gluon distribution has been adopted from 
Kwiecinski et al. ||l5[] . The parameters of the quarkonium 
bound state are the ones given by Eichten and Quigg . 
For the Xci state it is crucial that the gluons arc off-shell 
in fc I -factorization. 
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The color octet contribution 



Now we proceed with the calculation of J/'F produc- 
tion by Xc radiative decays adopting the colour octet 
mechanism. The infrared stability of higher order cor- 
rections to the cross section requires the existence of a 
color octet contribution, without fixing its size The 
Xc state can be written in a velocity expansion as [ p!o[ 

\Xe) - Oil) \qq ['P/] ) + 0(v) \qq Sf] g) + ■ ■ ■ . 

Following the formalism of the resulting cross sec- 

tion is then proportional to the color octet matrix el- 
ement (0| Og"^ (^S*!) |0) which has to be fitted to data. 
Using the results for the color singlet part and adding the 
color octet contribution we obtain as value of the color 
octet matrix element (0| Of ' (^Si) |0) (9.0 ± 2.0) x 
10~^. Comparing this with the result obtained in the 
coUinear factorization we find a suppression of the 

matrix element due to the flat Pj^-dependence of the color 
octet contribution by roughly one order of magnitude, re- 
sulting in a violation of the velocity scaling rules. These 
scaling rules are derived rigorously in the framework of 
non-relativistic QCD (NRQCD) It is, therefore, nat- 
ural to assume that the charm quark is simply not heavy 
enough for the velocity scaling rules of NRQCD to be 
valid. This is also suggested by other observations, see 
e.g. the very recent study Q. In contrast the descrip- 
tion of bottom systems in NRQCD should be more accu- 
rate. This shows the importance of a detailed analysis of 
bottomonia production in the /c^-factorization approach. 

Let us conclude. The /cj^-factorization approach re- 
lying on an unintegrated gluon distribution compatible 
with the small x behaviour of the structure function F2 
together with the BFKL NLLA fermion production ver- 
tices describes correctly Xc production in the central ra- 
pidity region. Whereas the standard coUinear factoriza- 
tion approach in NLO can describe the data in the TeV 
range only by introducing a dominant octet contribution. 



we have shown that in the fc_L-factorization approach such 
a contribution gives an improved description of the data 
but is suppressed by its P± behaviour. 

Our main conclusion is therefore that the correct way 
to improve the standard QCD calculations for quarko- 
nium production in the TeV range is to abandon the 
coUinear approximation. The contributions disregarded 
in the coUinear approximation of strong transverse mo- 
mentum ordering become essential in the small- a; range. 
The relative merits of the k± -factorization as the stan- 
dard approach for other processes in high energy hadronic 
collisions still has to be investigated. 
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